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Abstract 

Persistency is the property, for differential equations in M n , that solutions start- 
ing in the positive orthant do not approach the boundary. For chemical reactions 
and population models, this translates into the non-extinction property: provided 
that every species is present at the start of the reaction, no species will tend to be 
eliminated in the course of the reaction. This paper provides checkable conditions 
for persistence of chemical species in reaction networks, using concepts and tools 
from Petri net theory, and verifies these conditions on various systems which arise 
in the modeling of cell signaling pathways. 
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1 Introduction 



One of the main goals of molecular systems biology is the understanding of cell behavior 
and function at the level of chemical interactions, and, in particular, the characterization 
of qualitative features of dynamical behavior (convergence to steady states, periodic 
orbits, chaos, etc). A central question, thus, is that of understanding the long-time 
behavior of solutions. In mathematical terms, and using standard chemical kinetics 
modeling, this problem may be translated into the study of the set of possible limit 
points (the u> -limit set) of the solutions of a system of ordinary differential equations. 

Robustness 

A distinguishing feature of this study in the context of cell biology, in contrast to more 
established areas of applied mathematics and engineering, is the very large degree of 
uncertainty inherent in models of cellular biochemical networks. This uncertainty is due 
to environmental fluctuations, and variability among different cells of the same type, as 
well as, from a mathematical analysis perspective, the difficulty of measuring the relevant 
model parameters (kinetic constants, cooperativity indices, and many others) and thus 
the impossibility of obtaining a precise model. Thus, it is imperative to develop tools 
that are "robust" in the sense of being able to provide useful conclusions based only 
upon information regarding the qualitative features of the network, and not the precise 
values of parameters or even the forms of reactions. Of course, this goal is often not 
unachievable, since dynamical behavior may be subject to phase transitions (bifurcation 
phenomena) which are critically dependent on parameter values. 

Nevertheless, and surprisingly, research by many, notably by Clarke ^U]; Horn and 
Jackson |22| 128] , Feinberg ^Sl E31 E] , and many others in the context of complex balanc- 
ing and deficiency theory, and by Hirsch and Smith [SSI QUI and many others including 
the present authors [2J E21 El E| in the context of monotone systems, has resulted in 
the identification of rich classes of chemical network structures for which such robust 
analysis is indeed possible. In this paper, we present yet another set of tools for the 
robust analysis of dynamical properties of biochemical networks, and apply our approach 
in particular to the analysis of persistence in chemical networks modeled by ordinary 
differential equations. Our approach to studying persistence is based on the formalism 
and basic concepts of the theory of Petri nets. Using these techniques, our main results 
provide conditions (some necessary, and some sufficient) to test persistence. We then 
apply these conditions to obtain fairly tight characterizations in non-trivial examples 
arising from the current molecular biology literature. 

Persistency 

Persistency is the property that, if every species is present at the start of the reaction, 
no species will tend to be eliminated in the course of the reaction. Mathematically, this 
property can be equivalently expressed as the requirement that the cj-limit set of any 
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trajectory which starts in the interior of the positive orthant (all concentrations positive) 
does not intersect the boundary of the positive orthant (more precise definitions are 
given below). Persistency can be interpreted as non-extinction: if the concentration 
of a species would approach zero in the continuous differential equation model, this 
means, in practical terms, that it would completely disappear in finite time, since the 
true system is discrete and stochastic. Thus, one of the most basic questions that one 
may ask about a chemical reaction is if persistency holds for that network. Also from a 
purely mathematical perspective persistency is very important, because it may be used in 
conjunction with other tools in order to guarantee convergence of solutions to equilibria. 
For example, if a strictly decreasing Lyapunov function exists on the interior of the 
positive orthant (see e.g. [221 123 EH EH 113 Ell f° r classes of networks where this can be 
guaranteed), persistency allows such a conclusion. 

An obvious example of a non-persistent chemical reaction is a simple irreversible 
conversion A —>■ B of a species A into a species B; in this example, the chemical A 
empties out, that is, its time- dependent concentration approaches zero as t — > oo. This is 
obvious, but for complex networks determining persistency, or lack thereof, is, in general, 
an extremely difficult mathematical problem. In fact, the study of persistence is a classical 
one in the (mathematically) related field of population biology, where species correspond 
to individuals of different types instead of chemical units; see for example [T71 [7] and 
much other foundational work by Waltman. (To be precise, what we call "persistence" 
coincides with the usage in the above references, and is also sometimes called "strong 
persistence," at least when all solutions are bounded, a condition that we will assume in 
most of our main results, and which is automatically satisfied in most examples. Also, 
we note that a stronger notion, "uniform" persistence, is used to describe the situation 
where all solutions are eventually bounded away from the boundary, uniformly on initial 
conditions, see [HI EE]-) Most dynamical systems work on persistence imposes conditions 
ruling out phenomena such as heteroclinic cycles on the boundary of the positive orthant, 
and requiring that the unstable manifolds of boundary equilibria should intersect the 
interior, and more generally studying the chain-recurrence structure of attractors, see 

e.g. mi- 

Petri nets 

Petri nets, also called place/transition nets, were introduced by Carl Adam Petri in 
1962 j3JJ, and they constitute a popular mathematical and graphical modeling tool used 
for concurrent systems modeling [HOI EH]. Our modeling of chemical reaction networks 
using Petri net formalism is not in itself a new idea: there have been many works, at 
least since j32], which have dealt with biochemical applications of Petri nets, in particular 
in the context of metabolic pathways, see e.g. [TBI 123 123 123 121]; and especially the 
excellent exposition However, there does not appear to have been any previous 

work using Petri nets for a nontrivial study of dynamics. In this paper, although we do 
not use any results from Petri net theory, we employ several concepts (siphons, locking 
sets, etc.), borrowed from that formalism and introduced as needed, in order to formulate 
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new, powerful, and verifiable conditions for persistence and related dynamical properties. 



Application to a common motif in systems biology 

In molecular systems biology research, certain "motifs" or subsystems appear repeatedly, 
and have been the subject of much recent research. One of the most common ones is that 
in which a substrate So is ultimately converted into a product P, in an "activation" reac- 
tion triggered or facilitated by an enzyme E, and, conversely, P is transformed back (or 
"deactivated") into the original So, helped on by the action of a second enzyme F. This 
type of reaction is sometimes called a "futile cycle" and it takes place in signaling trans- 
duction cascades, bacterial two-component systems, and a plethora of other processes. 
The transformations of So into P and vice versa can take many forms, depending on how 
many elementary steps (typically phosphorylations, methylations, or additions of other 
elementary chemical groups) are involved, and in what order they take place. Figure ^ 
shows two examples, (a) one in which a single step takes place changing So into P = Si, 
and (b) one in which two sequential steps are needed to transform So into P = S2, with 
an intermediate transformation into a substance Si. A chemical reaction model for such 




Figure 1: (a) One-step and (b) two-step transformations 



a set of transformations incorporates intermediate species, compounds corresponding to 
the binding of the enzyme and substrate. (In "quasi-steady state" approximations, a 
singular perturbation approach is used in order to eliminate the intermediates. These 
approximations are much easier to study, see e.g. 2J.) Thus, one model for (a) would be 
through the following reaction network: 

E + S «-> ES -> £ + Si m 
F + Si <-> FSi -> F + S [ ' 

(double arrows indicate reversible reactions) and a model for (b) would be: 

E + S <-> ES -> F + Sx <-> ESi -> E + S 2 

F + S 2 «-> FS 2 -> F + Si <-> FSx -> F + S [ ' 

where "FSo" represents the complex consisting of E bound to So and so forth. 

As a concrete example, case (b) may represent a reaction in which the enzyme E 
reversibly adds a phosphate group to a certain specific amino acid in the protein So, 
resulting in a single-phosphorylated form Si; in turn, E can then bind to Si so as to pro- 
duce a double-phosphorylated form S 2 , when a second amino acid site is phosphorylated. 
A different enzyme reverses the process. (Variants in which the individual phosphory- 
lations can occur in different orders are also possible; we discuss several models below.) 
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This is, in fact, one of the mechanisms believed to underlie signaling by MAPK cascades. 
Mitogen- activated protein kinase (MAPK) cascades constitute a motif that is ubiquitous 
in signal transduction processes j2U EHl EZj in eukaryotes from yeast to humans, and 
represents a critical component of pathways involved in cell apoptosis, differentiation, 
proliferation, and other processes. These pathways involve chains of reactions, activated 
by extracellular stimuli such as growth factors or hormones, and resulting in gene expres- 
sion or other cellular responses. In MAPK cascades, several steps as in (b) are arranged 
in a cascade, with the "active" form 5*2 serving as an enzyme for the next stage. 

Single-step reactions as in (a) can be shown to have the property that all solutions 
starting in the interior of the positive orthant globally converge to a unique (subject to 
stoichiometry constraints) steady state, see |3j, and, in fact, can be modeled by monotone 
systems after elimination of the variables E and F, cf. p. The study of (b) is much 
harder, as multiple equilibria can appear, see e.g. [23111]. We will show how our results 
can be applied to test consistency of this model, as well as several variants. 

Organization of paper 

The remainder of paper is organized as follows. Section |21 sets up the basic terminol- 
ogy and definitions regarding chemical networks, as well as the notion of persistence, 
Section |3] shows how to associate a Petri net to a chemical network, Sections 0] and 
provide, respectively, necessary and sufficient conditions for general chemical networks, 
In Section ® we show how our results apply to the enzymatic mechanisms described 
above. 

2 Chemical Networks 

A chemical reaction network ("CRN", for short) is a set of chemical reactions TZi, where 
the index i takes values in TZ := {1, 2, . . . , n r }. We next define precisely what one means 
by reactions, and the differential equation associated to a CRN, using the formalism from 
chemical networks theory. 

Let us consider a set of chemical species Sj, j G {1,2, ...n s } := S which are the 
compounds taking part in the reactions. Chemical reactions are denoted as follows: 



where the cty and are nonnegative integers called the stoichiometry coefficients. The 
compounds on the left-hand side are usually referred to as the reactants, and the ones 
on the right-hand side are called the products, of the reaction. Informally speaking, the 
forward arrow means that the transformation of reactants into products only happens in 
the direction of the arrow. If also the converse transformation occurs, then, the reaction 
is reversible and we need to also list its inverse in the chemical reaction network as a 
separate reaction. 




(3) 
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It is convenient to arrange the stoichiometry coefficients into an n s x n r matrix, called 
the stoichiometry matrix T, defined as follows: 

[r]jr'j = Pij ~ (4) 

for alH E 1Z and all j E S (notice the reversal of indices). This will be later used in 
order to write down the differential equation associated to the chemical reaction network. 
Notice that we allow T to have columns which differ only by their sign; this happens when 
there are reversible reactions in the network. 

We discuss now how the speed of reactions is affected by the concentrations of the 
different species. Each chemical reaction takes place continuously in time with its own 
rate which is assumed to be only a function of the concentration of the species taking 
part in it. In order to make this more precise, we define the vector S — [Si, S 2 , ■ ■ ■ S ns ]' 
of species concentrations and, as a function of it, the vector of reaction rates 

R(S) := [Ri(S),R 2 (S),...R nr (S)}'. 

Each reaction rate Ri is a real-analytic function defined on an open set which contains 
the non-negative orthant + = M> s of M" -3 , and we assume that each Ri depends only on 
its respective reactants. (Imposing real-analyticity, that is to say, that the function Ri 
can be locally expanded into a convergent power series around each point in its domain, 
is a very mild assumption, verified in basically all applications in chemistry, and it allows 
stronger statements to be made.) Furthermore, we assume that each Ri satisfies the 
following monotonicity conditions: 

dRj(S) \ >0 ife^>0 
dSj \ = if aij = 0. [b) 

We also assume that, whenever the concentration of any of the reactants of a given 
reaction is 0, then, the corresponding reaction does not take place, meaning that the 
reaction rate is 0. In other words, if S^, ... , Si N are the reactants of reaction j, then we 
ask that 

Rj(S) = for all S such that [S h , . . . , S iN ] E dO + , 

where dO + = <91R> is the boundary of + in W N . Conversely, we assume that reactions 
take place if reactants are available, that is: 

Rj(S) > whenever S is such that [S h , . . . , S iN ] E int[R> ] . 

A special case of reactions is as follows. One says that a chemical reaction network 
is equipped with mass-action kinetics if 

Ri(S) = fci II S T° for alH = 1, . . . , n r . 

3=1 
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This is a commonly used form for the functions Ri(s) and amounts to asking that the 
reaction rate of each reaction is proportional to the concentration of each of its partic- 
ipating reactants. The results in this paper do not require this assumption; in a paper 
in preparation will we specialize and tighten our results when applied to systems with 
mass-action kinetics. 

With the above notations, the chemical reaction network is described by the following 
system of differential equations: 

S = TR(S). (6) 

with S evolving in + and where V is the stoichiometry matrix. 

There are several additional notions useful when analyzing CRN's. One of them is 
the notion of a complex. We associate to the network (jSJ) a set of complexes, Cj's, with 
i G {1,2,..., n c }. Each complex is an integer combination of species, specifically of the 
species appearing either as products or reactants of the reactions in (J3J). We introduce 
the following matrix T as follows: 

a u a 2 i ... a nr i (3u (3 2 i ... (3 nr i 
a 12 a 22 ■ ■ ■ «n r 2 /?i2 P22 ... (3 nr 2 

Then, a matrix representing the complexes as columns can be obtained by deleting from 
T repeated columns, leaving just one instance of each; we denote by T c G R™ sXnc the 
matrix which is thus constructed. Each of the columns of T c is then associated with a 
complex of the network. We may now associate to each chemical reaction network, a 
directed graph (which we call the C-graph), whose nodes are the complexes and whose 
edges are associated to the reactions (j2J). An edge (Cj, Cj) is in the C-graph if and only if 
Cj — > Cj is a reaction of the network. Note that the C-graph need not be connected (the 
C-graph is connected if for any pair of distinct nodes in the graph there is an undirected 
path linking the nodes), and lack of connectivity cannot be avoided in the analysis. (This 
is in contrast with many other graphs in chemical reaction theory, which can be assumed 
to be connected without loss of generality.) In general, the C-graph will have several 
connected components (equivalence classes under the equivalence relation "is linked by 
an undirected path to", defined on the set of nodes of the graph). 

Let I be the incidence matrix of the C-graph, namely the matrix whose columns are 
in one-to-one correspondence with the edges (reactions) of the graph and whose rows 
are in one-to-one correspondence with the nodes (complexes). Each column contains a 
— 1 in the i-th entry and a +1 in the j-th entry (and zeroes in all remaining entries) 
whenever (Cj, Cj) is an edge of the C-graph (equivalently, when Cj — > Cj is a reaction of 
the network). With this notations, we have the following formula, to be used later: 

r = r c x. (7) 

We denote solutions of © as follows: S(t) = ip(t,So), where So G + is the initial 
concentration of chemical species. As usual in the study of the qualitative behavior 
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of dynamical systems, we will make use of cu-limit sets, which capture the long-term 
behavior of a system and are defined as follows: 

co(Sq) := {S G + : <f(t n , So) — > S for some t n — > +00} (8) 

(implicitly, when talking about uj(S ), we assume that ip(t,S ) is defined for all t > 
for the initial condition So). We will be interested in asking whether or not a chemical 
reaction network admits solutions in which one or more of the chemical compounds 
become arbitrarily small. The following definition, borrowed from the ecology literature, 
captures this intuitive idea. 

Definition 2.1 A chemical reaction network (jUJ) is persistent if oo(So) H dO + = for 
each 5*0 G int(C + ). □ 

We will derive conditions for persistence of general chemical reaction networks. Our 
conditions will be formulated in the language of Petri nets; these are discrete-event 
systems equipped with an algebraic structure that reflects the list of chemical reactions 
present in the network being studied, and are defined as follows. 



3 Petri Nets 

We associate to a CRN a bipartite directed graph (i.e., a directed graph with two types 
of nodes) with weighted edges, called the species-reaction Petri net, or SR-net for short. 
Mathematically, this is a quadruple 

(V S ,V R ,E,W), 

where Vs is a finite set of nodes each one associated to a species, Vr is a finite set of 
nodes (disjoint from Vs), each one corresponding to a reaction, and E is a set of edges 
as described below. (We often write S or Vs interchangeably, or R instead of Vr, by 
identifying species or reactions with their respective indices; the context should make the 
meaning clear.) The set of all nodes is also denoted by V = Vr U Vs. 

The edge set E C V x V is defined as follows. Whenever a certain reaction Ri belongs 
to the CRN: 

jes jes 

we draw an edge from Sj G Vs to Ri G Vr for all S/s such that a^- > 0. That is, 
(Sj,Rj) G E iff aij > 0, and we say in this case that Ri is an output reaction for Sj. 
Similarly, we draw an edge from Ri G Vr to every Sj G Vs such that /3y > 0. That is, 
(Ri, Sj) G E whenever > 0, and we say in this case that Ri is an input reaction for 
Sj. 

Notice that edges only connect species to reactions and vice versa, but never connect 
two species or two reactions. 
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The last element to fully define the Petri net is the function W : E — > N, which 
associates to each edge a positive integer according to the rule: 

W(S J} Ri) = an and W(R U Sj) = . 

Several other definitions which are commonly used in the Petri net literature will be 
of interest in the following. We say that a row or column vector v is non-negative, and 
we denote it by v >z if it is so entry-wise. We write v y if v >z and v ^ 0. A 
stronger notion is instead u>0, which indicates f « > for all i. 

Definition 3.1 A P-semiflow is any row vector c y such that cT = 0. Its support is 
the set of indices {i G Vg : Cj > 0}. A Petri net is said to be conservative if there exists 
a P-semiflow c ^> 0. □ 

Notice that P-semiflows for the system © correspond to non-negative linear first 
integrals, that is, linear functions S i— > cS such that (d/dt)cS(t) = along all solutions 
of © (assuming that the span of the image of R(S) is M. nr ). In particular, a Petri net 
is conservative if and only if there is a positive linear conserved quantity for the system. 
(Petri net theory views Petri nets as "token-passing" systems, and, in that context, P- 
semiflows, also called place-invariants, amount to conservation relations for the "place 
markings" of the network, that show how many tokens there are in each "place," the 
nodes associated to species in SR-nets. We do not make use of this interpretation in this 
paper.) 

Definition 3.2 A T-semiflow is any column vector v y such that Tv = 0. A Petri 
net is said to be consistent if there exists a T-semiflow v ^> 0. □ 

The notion of T-semiflow corresponds to the existence of a collection of positive 
reaction rates which do not produce any variation in the concentrations of the species. 
In other words, v can be viewed as a set of fluxes that is in equilibrium (38.). (In Petri 
net theory, the terminology is "T- invariant," and the fluxes are flows of tokens.) 

A chemical reaction network is said to be reversible if each chemical reaction has an 
inverse reaction which is also part of the network. Biochemical models are most often 
non- reversible. For this reason, a far milder notion was introduced [221 EH EH 
A chemical reaction network is said to be weakly reversible if each connected component 
of the C-graph is strongly connected (meaning that there is a directed path between 
any pair of nodes in each connected component). In algebraic terms, weak reversibility 
amounts to existence of v ^> such that Xv = (see Corollary 4.2 of 16 ), so that in 
particular, using (JJJ), also Tv = T c 2v = 0. Hence a chemical reaction network that is 
weakly reversible has a consistent associated Petri net. 

A few more definitions are needed in order to state our main results. 
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Definition 3.3 A nonempty set E C Vs is called a siphon if each input reaction asso- 
ciated to E is also an output reaction associated to E. A siphon is a deadlock if its set 
of output reactions is all of Vr. A deadlock is minimal if it does not contain (strictly) 
any other deadlocks. A pair of distinct deadlocks Si and E2 is said to be nested if either 
Ei C E2 or E2 C Ei. □ 

For later use we associate a particular set to a siphon E as follows: 

L 2 = {x e + I Xi = i e E}. 

It is also useful to introduce a binary relation "reacts to", which we denote by >—>, and 
we define as follows: Si >— > Sj whenever there exists a chemical reaction IZf., so that 

les les 

with aki > 0, (3kj > 0. If the reaction number is important, we also write 

i ' y Oj 

(where k G 1Z). With this notation, the notion of siphon can be rephrased as follows: 
Z C S is a siphon for a chemical reaction network if for every S G Z and k G 1Z such 
that S k :={T eS:T AS} ^ 0, it holds S k n Z ^ 0. 



4 Necessary conditions 

Our first result will relate persistence of a chemical reaction network to consistency of 
the associated Petri net. 

Theorem 1 Let (EJ) be the equation describing the time- evolution of a conservative and 
persistent chemical reaction network. Then, the associated Petri net is consistent. 

Proof. Let So G int(C?+) be any initial condition. By conservativity, solutions satisfy 
cS(t) = cSq, and hence remain bounded, and therefore u(Sq) is a nonempty compact set. 
Moreover, by persistence, u)(Sq) H dO + = 0, so that R(S ) ^> 0, for all S G lj(Sq). In 
particular, by compactness of uj(Sq) and continuity of R, there exists a positive vector 
v ^> 0, so that 

R(Sq) >z v for all Sq G uj(Sq) . 

Take any So G u)(So). By invariance of u(So), we have R(<p(t,So)) ^ v for all {el 
Consequently, taking asymptotic time averages, we obtain: 

0= lim ^ ^ " ^° = hm I [ T TR( V (t,S ))dt (10) 

i ^ + 00 J i ^ + 00 J 
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(the left-hand limit is zero because <f(T, So) is bounded). However, 



1 



■r 



T 



R(<p(t, S )) dt y v 



for all T > 0. Therefore, taking any subsequence T n — > +00 so that there is a finite limit: 



lim — / R(ip(t, S )) dt = v y v . 

n->+oc T n J 

We obtain, by virtue of (|10|). that T w = 0. This completes the proof of consistency, since 



5 Sufficient Conditions 

In this present Section, we derive sufficient conditions for insuring persistence of a chem- 
ical reaction network on the basis of Petri net properties. 

Theorem 2 Consider a chemical reaction network satisfying the following assumptions: 

1. its associated Petri net is conservative; 

2. each siphon contains the support of a P-semiflow. 

Then, the network is persistent. 

We first prove a number of technical results. The following general fact about differ- 
ential equations will be useful. 

For each real number p, let signp := 1, 0, — 1 if p > 0, p — 0, or p < respectively, 
and for each vector x = (xx, . . . ,x n ), let sign a; := (signxi, . . . ,signx n )'. When x belongs 
to the closed positive orthant R", sign a; G {0, l} n . 

Lemma 5.1 Let / be a real- analytic vector field defined on some open neighborhood of 
IR™, and suppose that M™ is forward invariant for the flow of /. Consider any solution 
x(t) of x = f(x), evolving in R™ and defined on some open interval J. Then, signx(t) is 
constant on J. 

Proof. Pick such a solution, and define 




v > 0. 



Z := {i\xi{t) = for all t G J} . 
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Relabeling variables if necessary, we assume without loss of generality that Z = {r + 
1, . . . , to}, with < r < to, and we write equations in the following block form: 

y = g(v,z) 

z = h(y,z) 

where x' = (y', z')' and y(t) G W, z(t) G R n_r . (The extreme cases r = and r = to 
correspond to x = z and x = y respectively.) In particular, we write x! = (y', z')' for the 
trajectory of interest. By construction, z = 0, and the sets 

B l := {t\yi(t) = Q} 

are proper subsets of J, for each i G {1, . . . ,r}. Since the vector field is real-analytic, 
each coordinate function y^ is real-analytic (see e.g. [US], Proposition C.3.12), so, by the 
principle of analytic continuation, each Bi is a discrete set. It follows that 

r 

G:= J\(jA 

i=i 

is an (open) dense set, and for each t G G, y(t) G inter W + , the interior of the positive 
orthant. 

We now consider the following system on W: 

y = g(y, o) . 

This is again a real- analytic system, and W + is forward invariant. To prove this last 
assertion, note that forward invariance of the closed positive orthant is equivalent to the 
following property: 

for any y G W + and any i G {1, . . . , r} such that = 0, gi{y, 0) > 0. 

Since Wf is forward invariant for the original system, we know, by the same property 
applied to that system, that for any (y, z) G M™ and any % G {1, . . . , r} such that yi = 0, 
9i(y,z) > 0. Thus, the required property holds (case z = 0). In particular, inter R!j_ 
is also forward invariant (see e.g. |2j, Lemma III. 6). By construction, y is a solution of 
y — giVifyi y(t) G inter R!j_ for each t G G, Since G is dense and inter IR^_ is forward 
invariant, it follows that y(t) G inter W + for all t G J. Therefore, 

signx(t) = (l r , n _ r )' for all t G J 

where l r is a vector of r l's and n _ r is a vector of to — r 0's. ■ 

We then have an immediate corollary: 

Lemma 5.2 Suppose that Q C + is a closed set, invariant for (JHJ). Suppose that flHLz 
is non-empty, for some Z G S. Then, Q H Lz is also invariant with respect to ©• 
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Proof. Pick any Sq G QC\Lz- By invariance of Q, the solution ip(t, So) belongs to Q for all 
t in its open domain of definition J, so, in particular (this is the key fact), ip(t, S ) G + 
for all t (negative as well as positive). Therefore, it also belongs to Lz, since its sign is 
constant by Lemma f5. 11 I 

In what follows, we will make use of the Bouligand tangent cone TC^(K) of a set 
K C + at a point £ G + , defined as follows: 

TC^K) = ER n :3k n E K,k n ^^nd\ n \0: J-(fc n -f)^*;J . 

Bouligand cones provide a simple criterion to check forward invariance of closed sets (see 
e.g. [H|): a closed set K is forward invariant for (0) if and only if TR(£) G TC^(K) for all 
£ E K. However, below we consider a condition involving tangent cones to the sets Lz, 
which are not closed. Note that, for all index sets Z and all points £ in Lz, 

TQ (L z ) = {v G R n : v t = Vi G Z} . 

Lemma 5.3 Let Z C <S be non-empty and £ G be such that Ti?(£) G TC^(L Z ). 
Then Z is a siphon. 

Proof. By assumption Ti?(£) G TC^(L Z ) for some £ G L^. This implies that [T\R(£)]j = 
for all i E Z. Since = for all i G Z, all reactions in which Si is involved as a reactant 
are shut off at £; hence, the only possibility for [r j = is that all reactions in which 
Si is involved as a product are also shut-off. Hence, for all k G TZ, and all I G S so that 
Si >^Si, we necessarily have that i4(0 = 0. 

Hence, for all k G 7£ so that = {I G 5 : ^ >^>Si} is non-empty, there must exist an 
/ G Sk so that £z = 0. But then necessarily, / G Z, showing that Z is indeed a siphon. ■ 

The above Lemmas are instrumental to proving the following Proposition: 

Proposition 5.4 Let £ G + be such that w(£) D Lz ^ for some Z C <S. Then Z is a 
siphon. 

Proof. Let f2 be the closed and invariant set o>(£). Thus, by Lemma [5 .2^ the non-empty 
set Lz fl is also invariant. Notice that 

d[L z ] = \J L w . 

WDZ 

Moreover, Lw H Q is invariant for all W C S such that Lw H f2 is non-empty. Hence, 

c\[L z ] n n = \J[L w nn] 

WDZ 
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is also invariant. By the characterization of invariance for closed sets in terms of Bouli- 
gand tangent cones, we know that, for any rj £ cl[Lx] H Q we have 

TR( V ) £ TC v (n nd(L z )) C TC v (c\(L z )). 

In particular, for r\ £ L z PI Q (which by assumption exists), TR(rj) £ TC V (L Z ) so that, 
by virtue of Lemma 15.31 we may conclude Z is a siphon. I 

Although at this point Proposition 15.41 would be enough to prove Theorem El it is 
useful to clarify the meaning of the concept of a "siphon" here. It hints at the fact, 
made precise in the Proposition below, that removing all the species of a siphon from the 
network (or equivalently setting their initial concentrations equal to 0) will prevent those 
species from being present at all future times. Hence, those species literally "lock" a part 
of the network and shut off all the reactions that are therein involved. In particular, once 
emptied a siphon will never be full again. This explains why a siphon is sometimes also 
called a "locking set" in the Petri net literature. A precise statement of the foregoing 
remarks is as follows. 

Proposition 5.5 Let Z C S be non-empty. Then Z is a siphon if and only if cl(L z ) is 
forward invariant for (JHJ). 

Proof. Sufficiency: Pick £ £ L z ^ 0. Then forward invariance of cl(L^) implies that 
ri?(£) £ TC^(c\(L z )) = TC^(Lz), where the last equality holds since £ £ L z . It follows 
from Lemma f5. 31 that Z is a siphon. 

Necessity: Pick £ £ cl(L z ). This implies that £j = for all i £ Z U Z', where Z' G S 
could be empty. By the characterization of forward invariance of closed sets in terms of 
tangent Bouligand cones, it suffices to show that [ri?(£)]j = for all i £ Z, and that 
[r#(£)]i > for all i £ Z' whenever Z' ^ 0. Now by ©, 

[ri?(£)]i = 5> fei i4(£) -J>^(£) = £ftii2fc(£)-0 > 0, (11) 

which already proves the result for z £ Z' . Notice that the second sum is zero because 
if an > 0, then species % is a reactant of reaction I, which implies that i?;(£) = since 
£i = 0. So we assume henceforth that i £ Z. We claim that the sum on the right side of 
(fTTj) is zero. This is obvious if the sum is void. If it is non-void, then each term which 
is such that (3^ > must be zero. Indeed, for each such term we have that Rk{0 = 
because Z is a siphon. This concludes the proof of Proposition 15.41 I 

Proof of Theorem [2] 

Let £ £ int(C + ) be arbitrary and let Q denote the corresponding cu-limit set fl = w(£). 
We claim that the intersection of Q and the boundary of + is empty. 



14 




Figure 2: Petri net associated to reactions (j2J. 



Indeed, suppose that the intersection is nonemty. Then, would intersect Lz, for 
some / Z C 5. In particular, by Proposition 15.41 Z would be a siphon. Then, 
by our second assumption, there exists a non-negative first integral cS, whose support 
is included in Z, so that necessarily cS(t n ,£) — > at least along a suitable sequence 
t n — ► +oo. However, cS(t,£) — c£ > for all t > 0, thus giving a contradiction. ■ 

6 Applications 

We now apply our results to obtain persistence results for variants of the reaction (b) 
shown in Figure [T] as well as for cascades of such reactions. 

6.1 Example 1 

We first study reaction J2J). Note that reversible reactions were denoted by a in 
order to avoid having to rewrite them twice. The Petri net associated to (J2J) is shown 
if Fig. |21 The network comprises nine distinct species, labeled So, Si, S 2 , E, F, ES , 
ES\, FS2, FS\. It can be verified that the Petri net in Fig. |2] is indeed consistent 
(so it satisfies the necessary condition). To see this, order the species and reactions by 
the obvious order obtained when reading (j2J from left to right and from top to bottom 
(e.g., Si is the fourth species and the reaction E + Si — > ES\ is the fourth reaction). 
The construction of the matrix V is now clear, and it can be verified that Tv = with 
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v = [211211211211]'. The network itself, however, is not weakly reversible, since 
neither of the two connected components of (J2J) is strongly connected. Computations 
show that there are three minimal siphons: 

{E, ES , ESi}, 

{F,FS 1 ,FS 2 }, 

and 

{So, Si, S 2 , ES , ESi, FS 2 , FSi}. 

Each one of them contains the support of a P-semiflow; in fact there are three independent 
conservation laws: 

E + ES + ESi = const i, 

F + FS 2 + FSi = const 2 , and 

S + S 1 + S 2 + ES + ESt + FS 2 + FS 1 = const 3 , 

whose supports coincide with the three mentioned siphons. Since the sum of these three 
conservation laws is also a conservation law, the network is conservative. Therefore, 
application of Theorem |2] guarantees that the network is indeed persistent. 

6.2 Example 2 

As remarked earlier, examples as the above one are often parts of cascades in which 
the product (in MAPK cascades, a doubly-phosphorilated species) S 2 in turn acts as an 
enzyme for the following stage. One model with two stages is as follows (writing S 2 as 
E* in order to emphasize its role as a kinase for the subsequent stage): 

E + S <-> ES 
F + E* <-> FS 2 
E* + S^ <-> ES% 
+ <-> FS* 

The overall reaction is shown in Fig. El Note - using the labeling of species and reaction 
as in the previous example - that Tv = with v = [v[ v[ v[ v^]' and v± = [2 1 1 2 1 1]', and 
hence the network is consistent. There are five minimal siphons for this network, namely: 

{E, ES , ESi}, 

{F,FS 2 ,FS 1 }, 

{F*, FSI FSI}, 

{So, S±, S 2 , ESq, es*, fs 2 , fs*}, 

and 

{S , Si, E*, ESo, ESi, FS 2 , FS±, ESq, ES{}. 
Each one of them is the support of a P-semiflow, and there are five conservation laws: 



E + S! 
F + S 1 
E* + S\ 
F* + St 



ESi 
FS! 
ESf 
FS$ 



E + E* 
F + S 
E* + S$ 
F* + St. 



(12) 
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Figure 3: Petri net associated to reactions (|12j). 



E + ES + ESi = const i, 
F + FS 2 + FS X = const 2 , 
F* + FS* + FS$ = const 3 , 

S* + S* + S* + ES* + ES* + FS* + FS* = const 4 , 
and 

S + S 1 + E* + ES + ES 1 + FS 2 + FS 1 + ES* + ES* = const 5 . 

As in the previous example, the network is conservative since the sum of these conserva- 
tion laws is also a conservation law. Therefore the overall network is persistent, by virtue 
of Theorem El 



6.3 Example 3 

An alternative mechanism for dual phosphorilation in MAPK cascades, considered in |27j . 
differs from the previous ones in that it becomes relevant in what order the two phos- 
phorylations occur. (These take place at two different sites, a threonine and a tyrosine 
residue). The corresponding network can be modeled as follows: 

M + E ME -> M y + E M y E -> M 2 + E 

M + E ME* -> M t + E M t £ -> M 2 + £ , , 

M 2 + F <-> M 2 F -> M^ + F <-> M„F -> M + F ( °' 

M 2 + F .U>F* -> M, • /•' .!///•' -> M + F. 

See Fig. 01 for the corresponding Petri net. This network is consistent. Indeed, Tv = 
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Figure 4: Petri net associated to the network (|13j). 



for the same v as in the previous example. Moreover it admits three siphons of minimal 
support: 

{ /;. ME, ME*, M y E, M t E}, 

{F, M y F, M t F, M 2 F, M 2 F*}, 

and 

{M, ME, ME*, M y , M t , M y E, M t E, M 2 , M 2 F, M,F\ M t F, M y F}. 

Each of them is also the support of a conservation law, respectively for M,E and F 
molecules. The sum of these conservation laws, is also a conservation law and therefore 
the network is conservative. Thus the Theorem |2] again applies and the network is 
persistent. 
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